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We report the observation of macroscopic harmonic states in an optically induced confinement in a 
highly photoexcited semiconductor microcavity at room temperature. The spatially photomodulated 
refractive index changes result in the visualization of harmonic states in a micrometer-scale optical 
potential at quantized energies up to 4 meV even in the weak-coupling plasma limit. We characterize 
the time evolution of the harmonic states directly from the consequent pulse radiation and identify 
sequential multiple ~10 ps pulse lasing with different emitting angles and frequencies. 


Optical trapping of atoms is essential to realize atomic 
Bose-Einstein condensates [1]. Similarly, the 2D spatial 
confinement of excitons or exciton-polaritons by inhomo¬ 
geneous strains [2] , natural defects and potential fluctua¬ 
tions [3-6], nanofabrication [7-13], and optical potentials 
[14-17] in a semiconductor heterostructure/microcavity 
is considered conducive for the formation and control 
of the condensates of these quasiparticles. These dy¬ 
namic condensates can form a meta-stable state in a 
finite-momentum state [8, 18] and multiple spatial modes 
[19]. The optical visualization of macroscopic interacting 
quantum states in solid-state systems has been demon¬ 
strated in condensates of excitons (bound electron-hole 
pairs) and kindred quasiparticles. Such light-matter hy¬ 
brid condensates are typically formed at cryogenic tem¬ 
peratures in a photoexcited density much below the Mott 
transition where constituent quasiparticles can be con¬ 
sidered boson-like. Optically defined potentials enabled 
by the effective repulsive interactions of polaritons result 
in further real-time manipulation of these light-matter 
fluids in a steady state [15-17]. On the other hand, re¬ 
laxation oscillations with ^10-50 ps period caused by the 
interplay between reservoir feeding and Bose stimulation 
have been reported in localized condensates in natural 
potential fluctuations [5]. 

In this study, we report the optical visualization of 
dynamic quantized states in a highly photoexcited mi¬ 
crocavity at room temperature. Photoexcitation creates 
mainly free pairs of electrons and holes in a Ill-V-based 
semiconductor quantum well (QW) as a result of ther¬ 
mal ionization at high temperature [20, 21]. In the 
highly photoexcited microcavity studied here, the non- 
linearly photomodulated refractive index results in a siz¬ 
able effective cavity resonance shift, which enabls opti¬ 
cally induced confinement. In an optical confinement ini¬ 
tiated by spatially modulated nonresonant ps pulse ex¬ 
citation, sequential multiple-pulse lasing commences at 
several quantized energy levels. The transverse optical 
modes in a spatially modulated refractive index can also 
be understood as an optically induced potential for a 
fictitious quasiparticle (see Supplemental Material), e.g., 
correlated e-h pairs (eehp) resulting from the effective 


coupling to the cavity light field. Multiple-pulse las¬ 
ing from discrete states in a harmonic confinement (po¬ 
tential) manifests as a result of the time- and energy- 
dependent competition between the gain and reservoir 
carrier relaxation. 

The Fabry-Perot microcavity sample consists of a A 
GaAs cavity layer containing three sets of three In- 
GaAs/GaAs QWs embedded within GaAs/AlAs dis¬ 
tributed Bragg reflectors (DBRs) (see also Supplemen¬ 
tal Material). The sample is nonresonantly excited by 
a 2-ps pulse pump laser at Ep = 1.58 eV at room tem¬ 
perature. The pump energy is about 250 meV above the 
QW band gap {Eg ^ 1.33 eV) and 170 meV above the 
cavity resonance {Ec ~ 1.40-1.41 eV). The high-density 
e-h plasmas of a density ^ 1 — bx 10^^ cm“^ per QW per 
pulse are formed momentarily after non-resonant pulse 
excitation as a result of rapid (<10 ps) energy dissipation 
through optical phonons. The radiative recombination 
rate of these e-h carriers in the reservoir is suppressed 
because the cavity resonance Ec is detuned to ^ 70-80 
meV above the QW bandgap E'g^ i.e., the elhhl tran¬ 
sition between the first quantized electron and (heavy- 
)hole states in a QW. Below the lasing threshold, the 
high-density e-h plasmas in the microcavity are subject 
to non-radiative loss with a long decay time (>500 ps) 
(Supplemental Fig. SI). Therefore, the chemical poten¬ 
tial of the e-h plasmas (/io) appears to be stationary 
within ^100 ps after pulse excitation. The bare cav¬ 
ity resonance Ec of the sample studied here is close to 
Eg, the e2hh2 transition between the second quantized 
electron and hole levels of InGaAs/GaAs QWs. 

When jiiQ advances toward Ec, the average refractive 
index near Ec can be considerably modified, and this 
results in a sizable blueshift of the effective cavity reso¬ 
nance (F') and consequent emission energy (Supplemen¬ 
tal Fig. S2a). In general, F' increases with the photoex¬ 
cited density and can be seen as an effective potential V 
for eehps (Supplemental Material section S3). Therefore, 
a quasi-stationary confining potential V (r) for eehps can 
be established by a ring-shaped spatial distribution of 
photoexcited carrier density. In our experiments, we use 
a double-hump-shaped beam to fixate the orientation of 
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FIG. 1. Visualization of the macroscopic harmonic states, (a) Intensity image of the ring-shaped pump laser beam, (b) 
Photoluminescence (PL) image under a pump flux of about 1.3 Pth^ where the threshold pump flux Pth = 1.8 x 10® photons per 
pulse. The white dashed line represents the intensity peak of the pump. PL emerges at the center with a minimal overlap with 
the annular pump laser beam, (c-d) r-space imaging spectra at P = 0.8 Pth and 1.3 Pth- The black dashed line represents the 
harmonic confining potential V{x), whereas the white lines represent the spatial probability distributions of the lowest three 
states of a corresponding harmonic oscillator, (e-f) k-space imaging spectra. The energy splitting is hw ^ 2 meV, consistent 
with the quantized energy of a quantum oscillator for a particle with mass m* = 3 x 10~^ me, as determined by the E vs. k\\ 
dispersion (doted grey line). The quantized modes spectrally blue shift about 1 meV from P = 0.8 to 1.3 Pth, whereas the 
quantized energy splitting remains the same. The potential and spectral shifts are due to a density-dependent increase in the 
chemical potential of the high-density e-h plasma in the reservoir. 


the optically defined potential (Fig. la-b). 

Real-space (r-space) imaging spectra provide direct vi¬ 
sualization of the optical potential. Fig. Ic-d shows the 
r-space imaging spectra of a narrow cross-sectional stripe 
across the trap. A nearly parabolic potential well of ^10 
meV across 3 /im is revealed. Such a quasi-ID harmonic 
potential, V{x) = 1/2 is spontaneously formed un¬ 
der 2 ps pulse excitation. The resultant radiation ap¬ 
pears in-between the two humps, a few micrometers away 
from the pump spot. The r-space intensity profiles agree 
with the probability distributions of the quantized states 
of a harmonic oscillator (harmonic states). The devia¬ 
tions of the actual potential from a perfectly harmonic 
trap result in slightly asymmetric luminescence inten¬ 
sity distributions. These standing wave patterns form 
in the self-induced harmonic optical potential when a 
macroscopic coherent state emerges. At a critical density, 
the occupation number of cehps {rii) in these quantized 
harmonic states approaches unity when the conversion 
from the reservoir {Neh) overcomes the decay of the har¬ 
monic states. Consequently, when a stimulated process 
(oc Neh'^i) prevails, rii and the resultant radiation in¬ 
crease nonlinearly by a few orders of magnitude with the 
increasing N^h- 

An even more regular pattern appears in k-space imag¬ 
ing spectra (Fig. le-f). The E vs. k\\ dispersion mea¬ 
sured below the threshold allows the direct measurement 
of the effective mass m* = 3 x lO^rUe, where rrie is the 
electron rest mass. Moreover, the strength of the opti¬ 
cally defined harmonic potential a can be tuned through 
a variation in the spatial distance between two humps 
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FIG. 2. Quantized states in optically controlled con¬ 
fining potentials. K-space imaging spectra under below- 
threshold (a-b) and above-threshold (d-e) for two double¬ 
hump-shaped pump beams with peak-to-peak distances of 5 
/im and 3 /im, respectively. For comparison, the k-space imag¬ 
ing spectra under a flat-top pump beam are shown in (c) and 
(f). The corresponding r-space images and spectra are shown 
in Supplemental Fig. SI. 


of the pump beam (Fig. 2). The energy quantization 
(huj) varies with while the emission patterns 

evolve into the probability distributions for a particle 
with an effective mass m* in V{x) (Fig. 2 and Supple¬ 
mental Fig. S3). 

Temporally, these harmonic states emerge sequentially 
and display distinct density-dependent dynamics. In Fig. 
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FIG. 3. Dynamics, (a) Time-dependent luminescence in k-space at P = 1.0, 1.1, 1.2 and l.QPth- The P 3 , P 2 and Ei states 
appear sequentially with the increasing pump flux. The rise times decrease with the increasing pump flux for all states, (b) 
Time-dependent spectra in r-space. The false color represents normalized intensities. 


3a, we study the time evolutions of three harmonic states 
in k-space. The corresponding time-integrated imaging 
spectra in k-space are shown in Supplemental Fig. S4. 
The energy relaxation of these three states is revealed in 
the time-resolved spectra in Fig. 3b. At a critical pho- 
toexcited density, the high-energy P3 state arises ^ 25 
ps after pulse excitation and lasts for ^ 20 ps. The cor¬ 
responding pump flux is defined as the threshold (Pth)- 
With the increasing pump flux, the E 2 state emerges ^ 
25 ps after E^ at 1.1 Pfh^ whereas the ground Ei state 
appear 50 ps after E 2 at 1.2 Pt^. In the optically induced 
harmonic confinement studied here, the effective cavity 
resonance E'^ decreases with time as a result of the de¬ 
cay of reservoir carriers. However, the conversion from 
the reservoir to specific confined Ei state can be efficient 
only when E'^ is resonant with Ei. Such a temporally 
decrease in E'^ results in a time- and energy-dependent 
conversion efficiency for these confined harmonic states 
(Fig. SI) and consequent multiple-pulse lasing above the 
critical density threshold. 

We further analyze the emission flux and energy of 
these harmonic states by using time-integrated spectra 
measured with increasing photoexcited density (Fig. 4 
and Supplemental Fig. S2). Far below the threshold 
(p < 0.4 Pth)i emission is dominated by luminescence 
from GaAs spacer layers. When the pump flux is in¬ 
creased, the emissions from the cehp states become in¬ 
creasingly dominant, and the E^ state eventually lases 
at the threshold. The emission fluxes of all three states 
increase nonlinearly by more than two orders of mag¬ 


nitude across a threshold and then reach a plateau at 
a saturation density (Fig. 4a). On the other hand, the 
emission energy of these three states increases to a con¬ 
stant with the increasing pump flux (Fig. 4b). The en¬ 
ergy spacing huj only increases slightly with density. The 
spectral linewidths increase slightly for the Ei state but 
by about a factor of 10 for the E^ state. Next, we study 
the density-dependent dynamics. Fig. 4c shows the rise 
times and pulse durations for and Ei. The prod¬ 

uct of the variances of the spectral linewidth (AE) and 
the pulse duration (At) is found to be close to that of a 
transform-limited pulse: > Ah and ~ l/i for the E^/Ei 
and E 2 states, respectively. These harmonic states are 
macroscopically coherent states with finite phase and in¬ 
tensity fluctuations induced by interactions. 

We use a rate-equation model to describe the dy¬ 
namic formation of quantized states in an optically 
defined harmonic conflnment (Supplemental Material 
Sec. S4Supplemental Material). This phenomenolog¬ 
ical model reproduces qualitatively the dynamics and 
integrated emission flux of the harmonic states when 
photoexcited density is varied (Supplemental Figs. S5, 
S6, and S7). Non-equilibrium polariton condensates 
have been modeled by a modified Gross-Pitaevskii (GP) 
(or complex Ginzberg-Landau [cGL]) equation that ac¬ 
counts for the finite lifetime of polaritons [22, 23]. How¬ 
ever, cGL-type equations are inapplicable for the multi¬ 
ple dynamic states examined in this study. Additionally, 
the formation of a BEC-like condensate that underpins 
the GP- or cGL-type equation is not necessarily justified 
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FIG. 4. Density dependence, (a) Temporally and spec¬ 
trally integrated emission flux vs. pump flux. All three modes 
display non-linear increases in intensity by more than two or¬ 
ders of magnitude, and saturate at 1.1, 1.2 and 1.3 Pth, respec¬ 
tively. (b) Peak energy (solid shapes) and linewidths 2AE 
(error bars) vs. pump flux. These states spectrally blue shift 
by 1 to 4 meV. The spectral linewidths (AE) and pulsewidths 
(At) are reciprocal with a product of AE x At ^ Ah (h) for 
Es and Ei {E 2 ), which is closed to the uncertainty (Fourier- 
transform) limit, (c) Rise time vs. pump flux for the three 
states El (blue), E 2 (red) and E 3 (black). The error bar 
represents 2At. 


in our room-temperature experiments. 

Transverse light-held patterns and conhned optical 
modes have been identihed in nonlinear optical sys¬ 
tems [24], vertical-cavity surface-emitting lasers (VC- 
SELs) [25, 26], and microscale photonic structures [27]. 
In principle, the multiple transverse mode lasing in a 
high-density e-/i-plasma described in the present study 
can be modeled by a self-consistent numerical analy¬ 
sis with Maxwell-Bloch equations developed for conven¬ 
tional semiconductor lasers [28, 29], provided that the 
strong optical nonlinearities induced by Coulomb many- 


body effects, such as screening, bandgap renormalization, 
and phase-space hlling, are all considered. For example, 
one can consider the formation of index-guided multiple 
transverse modes as a result of an optically-induced re¬ 
fractive index reduction (Anc(x)). The cavity resonance 
shift (SE) can be estimated from 6Ec/Ec = —Aric/nc, 
where ric is the effective refractive index averaged over 
the longitudinal cavity photon mode which spans over 
^ /im in the growth direction. Therefore, a cavity res¬ 
onance shift 5E ^10 meV (Figs. 1 and 2) corresponds 
to \Anc/nc\ r\j 1%. Such a signihcant refractive index 
is probable with resonance-enhance optical nonlinearity 
from the e-h correlation [30-33] or carrier-induced change 
in refractive index at a high carrier density (of ^ 10^^ 
cm“^ or more) [34-36]. To uncover the microscopic for¬ 
mation mechanisms of such a sizable spatial modulation 
in refractive index or equivalent effective harmonic con¬ 
finement under pulse excitation, further characterizing 
the photoexcited density distribution with the use of 
other ultrafast spectroscopic techniques, such as a pump- 
probe spectroscopy, is necessary. 

We identify sequential multiple ^10 ps pulse lasing in 
an optically induced harmonic confinement in a semicon¬ 
ductor microcavity at room temperature. Laser radia¬ 
tion emerges at the quantized states of an optically in¬ 
duced harmonic confinement. The lasing frequency, rise 
time, pulse width, and radiation angle can be controlled 
through a variation in the photoexcited density or opti¬ 
cal pump spot dimensions in real time. The sample has 
a composition structure similar to those used for stud¬ 
ies of polariton condensates/lasers and the widely used 
VCSELs. In the highly photoexcited microcavity stud¬ 
ied here, harmonic confinement can be optically induced 
and controlled as a result of photomodulated refractive 
index changes even in the weak-coupling plasma limit. 
Our demonstration of macroscopic harmonic states in 
a room-temperature optically induced harmonic confine¬ 
ment improves our understanding of nonlinear laser dy¬ 
namics and should stimulate studies on emergent ordered 
states near the Fermi edge of a high-density e-h plasma 
in a semiconductor cavity [37-39]. 
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S.l. MATERIALS AND METHODS 


Sample. The microcavity sample is grown on a semi-insulating (lOO)-GaAs substrate with the molecular beam 
epitaxy method. The top (bottom) DBR consists of 17 (20) pairs of GaAs(61-nm)/AlAs (78-nm) A/4 layers. The 
central cavity layer consists of three sets of three Ino. 15 Gao. 85 As/GaAs (6 nm /12 nm) quantum wells each, positioned 
at the anti-nodes of the cavity light field. The structure is undoped and contains a A GaAs cavity sandwiched by 
DBRs. The bare cavity resonance Ec ~ 1.40 eV at room temperature. The QW bandgap {Eg ^1.33 eV) is tuned 
through a rapid thermal annealing process (at 1010°G-1090°G for 5-10 s), in which the InGaAs QW bandgap blueshifts 
because of the diffusion of gallium ions into the MQW layers. With increasing photoexcited density, the chemical 
potential jiQ of the e-h plasma can increase up to ^80 meV above Eg. The cavity quality factor Q is about 4000-7000, 
corresponding to a cavity photon lifetime of ^2 ps. 


Optical excitation. The front surface of the sample is positioned at the focal plane of a high-numerical-aperture 
objective (N.A. = 0.42, 50x, effective focal length: 4 mm). The holographic beam shaping is achieved with a reflected 
Fourier transform imaging system that consists of a 2D phase-only spatial light modulator (SLM), a 3 x telescope, a 
Faraday rotator, a polarizing beam splitter, and the objective. The light fields at the SLM and the sample surface 
form a Fourier transform pair. The 2D SLM (1920 x 1080 pixels, pixel pitch = 8 /im) enables us to generate arbitrary 
pump geometries with a ~2 jim spatial resolution at the sample surface by using computer-generated phase patterns. 
The polarization properties of pump and luminescence are controlled/analyzed with liquid-crystal devices. The pump 
flux is varied by more than two orders of magnitude with a liquid-crystal-based attenuator. 


Thermal management. At a high pump flux, the steady-state incident power transmitted to the sample at the 76 
MHz repetition rate of the laser would exceed 50 mW, and significant thermal heating results. Thermal heating can 
inhibit laser action and lead to spectrally broad redshifted luminescence. To suppress thermal heating, we temporally 
modulate the 2 ps 76 MHz pump laser pulse train with a duty cycle (on/off ratio) < 0.5% by using a double-pass 
acousto-optic modulator system. The time-averaged power is limited to below 1 mW for all experiments. 

R-space and k-space imaging spectroscopy. We measure the angular, spectral, and temporal properties of 
luminescence in the reflection geometry. The angle-resolved (k-space) luminescence images and spectra are measured 
through a Fourier transform optical system. A removable / = 200 mm lens enables the projection of either the k-space 
or r-space luminescence on to the entrance plane or slit of the spectrometer. Luminescence is collected through the 
objective, separated from the reflected specular and scattered pump laser light with a notch filter, and then directed 
to an imaging spectrometer. A lasing mode with a spatial diameter ^8 /rm is isolated for measurements through a 
pinhole positioned at the conjugate image plane of the microcavity sample surface. The spectral resolution is ^0.1 
nm (150 yueV), which is determined by the dispersion of the grating (1200 grooves/mm) and the entrance slit width 
(100-200 jim). The spatial (angular) resolution is ^ 0.3 /im (6 mrad) per GGD pixel. 
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S.2. SUPPLEMENTARY FIGURES 



Fig. SI. Below-threshold dynamics, (a) Time-dependent spectrally and spatially integrated luminescence at P = 0.9 Pth- 
The luminescence decays with a > 500 ps time constant largely because of the nonradiative loss, (b) Spatial luminescence 
distribution integrated up to a delay of 50 ps after pulse excitation, revealing a double-hump-shaped profile, (c) Time- 
dependent luminescence spectra at P = 0.9 Pth- The false color represents the square root of the luminescence intensity, (d) 
Time-dependent luminescence intensity for the P 3 (black), E 2 (red) and Pi (blue) harmonic states. The intensity is integrated 
over a spectral range of 1 meV centering the energies, as indicated by the dashed lines in (c). 
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Fig. S2. Pump-flux-dependent resonance energy and polarized luminescence spectra, (a) Resonance energy (the 
energy of emission mode at /cy 0 in K-space spectral mapping) versus photoexcited density under a tightly-focused circular 
pumping beam (radius ~ 1 /xm). The spectral peak can be considered as the effective cavity resonance and displays a sizable 
spectral blueshift with increasing the photoexcited density. Such an energy shift with photoexcited density is the basis of 
optically induced harmonic confinement (potential) using a spatially modulated pump beam profile. D ^ 4 x 10^^ cm“^ per 
quantum well per pulse. (b,c) Normalized spectra integrated over k-space as a function of pump flux in the double-hump-shaped 
pumping case: Co-circularly (b) and cross-circularly (c) polarized component. Quantized harmonic states are formed above 
~ O.SPth- Below the threshold, the luminescence from all states is unpolarized. Above the threshold, the E 3 state is highly 
circularly polarized, where as the Ei state has a diminishing circular polarization. The fourth quantized state {E 4 ) is weakly 
confined. 
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Fig. S3. Quantized states vs. the trap size, (a-c) Optical images of pump beam profiles at the front sample surface for 
two double-hump-shaped and one flat-top focal spots. The false color represents the photoexcited density (cm“^) per quantum 
well per pulse at 0.8 Pth for (a) and at 0.4 Pth for (b). Pth = 1-8 x 10® per pulse, as dehned in the main text. The lasing 
thresholds are 0.8 Pth for (a) and 0.4 Pth for (b). (d-f) Real-space imaging spectra measured below the threshold. The black 
dashed line shows the effective harmonic potential, whereas the white lines represent the spatial probability distributions of the 
lowest few states of a corresponding quantum oscillator. The quantized energy splittings for (d) and (e) are 2.2 and 4.0 meV, 
respectively. The false color represents the square root of intensity, (g-h) Real-space imaging spectra under above-threshold 
pump, (i) The measured wave functions of the two lowest states E 2 and Ei at 0.4 Pth (solid blue lines) and 0.5 Pth (solid 
red lines) as shown in (e) and (h). The ideal single-particle wave functions of the corresponding quantum simple harmonic 
oscillator (SHO) are represented by the dashed lines. 
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Fig. S4. K-space imaging spectra. The energy vs. in-plane momentum k\\ dispersion corresponding to the data set in the 
main paper with the increasing pump flux from 0.9 to 1.3 Pth- The quantized E 3 , E 2 , and Ei states lase in sequence with the 
increasing pump flux. F^a-state lasing commences hrst, followed by F^ 2 -state lasing at about 1.2 Pth, and then F^i-state lasing 
appears at 1.3 Pth and dominates over the other two states at a higher pump flux. 
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S.3. OPTICALLY INDUCED POTENTIAL AND REFRACTIVE INDEX CHANGES 

The transverse optical modes confined by a spatially modulated refractive index can also be understood as an 
optically induced potential for a fictitious quasiparticle. 

Considering the electromagnetic fields in a planar Fabry-Perot cavity with a medium of a 2D spatially modulated 
effective refractive index profile n(r), we can express the Helmholtz equation for the electric field in the following 
simplified form: 


d^F{r,ip,z) 

dr‘^ 


+ v?{r) ■ kl ■ F{r, ip,z) =0 


( 1 ) 


The electric field F in the cavity can be decomposed into lateral and transverse components, and approximated as 
F(r, ip, z) oc '0(r) •exp[—in(r) kz{r) • z)], where = /c^(r) +/cy (r), and k\\ (r) and kz(r) are transverse and longitudinal 
wavenumbers in the vacuum, respectively. For simplicity, we consider the ID case and reduce Eq. 1 to: 


d^F(x,z) 

dx^ 


+ n^{x) • [kl - k‘l{x)] • F{x, z) = 0 


( 2 ) 


Applying the effective-mass approximation for small k\\, we express the in-plane dispersion Ec{k\\) in terms of an 


771 (nhkz) 

■^11 = 


(n hk \\) 


where m* = 


hn^kz 


effective mass m* for the excitation: Ec{k\\) = E^o i — 2 m* ' 2m* ’ uv _ ^ 

Eq. 2 can then be mapped to a time-independent Schrdinger equation by replacing the electric field amplitude 
with the probability wavefunction 0 of a hypothetical quasi-particle with an effective mass m* in the presence of a 
potential V{x), thus: 


d‘^(j){x) 


+ [V{x) - E,] . (t){x) = 0 


2m* dx‘^ 

Replacing the effective mass for an optical cavity mode with longitudinal wavenumber kz, we obtain 

d‘^(j){x) in?kz 


dx^ 


he 


[E, - V{x)] . 0(x) = 0 


Comparing Eqs. 2 and 4 gives 


( 3 ) 

( 4 ) 


V{x) = he - kz{x) = Eco{x) 


( 5 ) 


Therefore, Eco{x) and E\\ correspond to the potential V{x) and the kinetic energy of a quasi-particle of mass m*. 
Considering Eco{x) oc where Lc is the effective cavity length, we can relate the potential to the spatially 

[1 + « E^o (1 - where h and V 


n(x) Lc 

modulated refractive index as follows: V{x) = V • 

are the refractive index and potential at x = 0. Thus, the spatial modulation of the potential {AV{x)) and refractive 
index {An{x)) have a one-to-one correspondence: 


AV{x) 
V 


An{x) 


S.4. THEORETICAL MODELING 

As discussed in section S.3, the spatial modulation of the refractive index can be considered a spatially dependent 
effective potential for a quasiparticle with an effective mass associated with the bare cavity resonance and mean 
refractive index. The multiple transverse optical modes are equivalent to harmonic states in harmonic confinement 
(potential). In this section, we use a set of rate equations that consider the formation of macroscopic harmonic states 
via the stimulation of a fictitious quasiparticle (correlated e-h pairs [cehps]) in the presence of an effective harmonic 
potential. We consider the temporal evolution of the reservoir distribution NR{x,t) and quantized sates ni{t). The 
e-h pairs photoexcited nonresonantly by a 2 ps pulse laser cool down rapidly (<5 ps) to the band edge. The cooled 
carriers in the reservoir Nji{x,t) are subject to nonradiative loss {Tnr) and thus result in a slow decrease (^0.1 
meV/ps) in the chemical potential po of the reservoir carriers. In this model, the bare cavity resonance Ec is fixed. 
By contrast, the effective cavity resonance E'^{t) blueshifts with the increasing photoexcited density (Fig. S2) and 
redshifts with time as a result of the temporal decay of reservoir carriers (Fig. S5). Standing waves (macroscopic 
coherent states of cehps) are formed in an optically induced harmonic confinement (potential V{x) oc x‘^) associated 
with photomodulated refractive index changes initiated by a double-humped pump beam profile. Such a harmonic 
potential is quasi-stationary within ^100 ps after the pulse excitation as a result of ^ns decay of the reservoir carriers. 
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A fraction of the reservoir carriers [Neh{x,t) = /3Nji{x,t)] near the effective cavity resonance E'^ are considered as 
cehps that can be scattered into the harmonic states via a spontaneous or stimulated process. Those cehps with 
energy nearly in resonance with a specific Ei state can be spontaneously scattered into the Ei state at a rate Wg. 
Above the threshold, the population of Ei states increases nonlinearly as a result of stimulation (oc Wss Eeh The 
consequent leakage photons from the decay of these confined cehps at a rate associated with the cavity photon decay 
rate Tc are then measured experimentally. The dynamics of the quantum oscillator of cehps can then be described by 
the following set of coupled rate equations: 


dNR{x,t) 

dt 

dni{t) 

dt 


Gp(x, t)P -Tnr Nr(x, Q Neh{x, t) [Wgg ni{t) + Ws] , 

i=l 

j dx7]i(iNeh{x,t) [Wss ni{t) + Ws] - Tc ni{t). 


(6) 

( 7 ) 


The spatial generation rate G{x,t) follows the 2 ps Gaussian temporal and double-hump-shaped beam profile 
(Fig. S5c and Fig. S6c). P is the pump fiux. Q is a scaling factor associated with (3 for a specific Ei state. r]i = 
exp \—{E'^ — Ei)‘^/{2 is the coupling efficiency of the reservoir cehp^ to the Ei state, where /S.Ei is the spectral 

linewidth of the Ei state. 


The fitting parameters are as follows: (3 = 0.0125, F^r = 1/700 ps Fc = 0.5 ps IFs = 1.0 x 10 ^ ps Wss 
= 1.0 X 10-2 ps-\ Cl = 0.55, C 2 = 0.67, Cs = 0.80, AEi = AE = 1 meV, Ei = 1.401 eV, E 2 = 1.403 eV, and Es = 
1.405 eV. In this simulation, we neglect the spin degree of freedom and only consider the density regime above 0.8 
and — Ei\ < 15 meV {Nj^ ^ 500, corresponding to the occupation number at P = Pth)- We approximate 

E'^{x^t) = PcO + 9 Neh{x^t)^ where PcO and g are determined by the following initial conditions: E'^{x = 0,t = 0) = 
1.4 eV and E'^{x = = 0) = 1.407 eV. 


This model reproduces the time evolution of the quantized states with increasing photoexcited density (Fig. S6). 
Below the threshold, the scattering of the reservoir cehps into the harmonic states occurs through a spontaneous decay 
process. The radiative quantum efficiency is low (< 10“^) because of the dominant non-radiative recombination loss. 
Above the threshold, the quantum efficiency increases by two to three orders of magnitude as a result of stimulation. 
This model reproduces qualitatively the rise times and the time-integrated emission fluxes of the multiple pulse lasing 
from the quantized states when the pump flux is varied (Fig. S7). 


We note a few shortcomings of this simplified model: (1) The model does not consider the local energy shift and 
the evolution of potentials induced by the interactions of the cehps because the density-dependent E'^{Nr, N^h^'^ii) is 
not well determined. As a result, the density-dependent energy shifts and linewidths cannot be reproduced. (2) The 
model does not consider carrier diffusion, which affects dynamics at a density far above the threshold. (3) Phase and 
intensity fluctuations are neglected. 




